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Abstract 

The single transverse spin asymmetries in heavy quark and anti-quark production from the 
quark-antiquark annihilation channel contribution is studied by taking into account the initial and 
final state interactions effects. Because of the different color charges, the final state interaction 
effects lead to about a factor of 3 difference in the spin asymmetry for heavy quark over that for 
the anti-quark in the valence region of low energy pp collisions. The experimental study of this 
model-independent prediction shall provide a crucial test for the underlying mechanism for the 
single spin asymmetry phenomena. 
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1. Introduction. Single-transverse spin asymmetry (SSA) is a novel and long stand- 
ing phenomena in hadronic reactions [US], and has attracted intensive interests from both 
experiment and theory sides in the last few years. It is a transverse spin dependent ob- 
servable, defined as the ratio of the cross section difference when we flip the transverse 
spin of one of the hadrons involved in the scattering over the sum of the cross sections, 
An oc (da(S±) — da(—S±))/(da(S±) + da(—S±)). The experimental observation of SSAs in 
semi-inclusive hadron production in deep inelastic scattering (SIDIS), in inclusive hadron 
production in pp scattering at collider energy at RHIC, and the relevant azimuthal asym- 
metric distribution of hadron production in e + e~ annihilation have motivated the theoretical 
developments in the last few years. It was also argued that the SSA_phenomena is closely 
related to the orbital motion of quarks and gluons in the nucleon 



Theoretically, it has been realized that the initial 



/final state interactions are crucial to 



a. 



These effects also invalidate the 



leading to a nonzero SSA in high energy scattering 
naive-time-reversal invariance argument for the non-existence of the so-called 

Sivers function JiJ, which describes a correlation between the intrinsic transverse momen- 
tum of parton and the transverse polarization vector of the nucleon. The initial/final state 
interaction effects introduce a pro cess- dependence in the SSA observables. For example, the 
SSA in the SIDIS process comes from final state interaction whereas that in the Drell-Yan 
lepton pair production comes from the initial state interaction. The consequence o 
difference is that there is a sign change between the SSAs for these two processes 
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It is of crucial to test this nontrivial QCD predictions by comparing the SSAs in these two 



processes. The Sivers sing 
HERMES collaboration 



e spin asymmetry in SIDIS process has been observed by the 
121, and the planned Drell-Yan measurement at RHIC and other 

n 

facility will test this prediction [131 ] . 

In this paper, we study another interesting probe for the initial/final state interaction 
effects: the SSAs in heavy quark and antiquark production in hadronic process. Because 
the heavy quark and antiquark can be detected by their decay products, their SSAs can be 
measured separately. The heavy quark and antiquark produced in short distance partonic 
processes will experience different final state interactions with the nucleon remanet due to 
their different color charges, and therefore the SSAs for heavy quark and antiquark will be 
different. Detailed calculations in the following show that the difference could be as large 
as a factor of 3 if the quark-antiquark channel contribution dominates. This is a direct 
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FIG. 1: Initial/final state interactions contributions to the SSAs in heavy quark and antiquark 
productions: (a) initial state interaction with color strength 1/2N C ; (b) final state interaction on 
the quark line with color strength (N^ — 2)/2iV c ; (c) final state interaction on the antiquark line 
with color strength 2/2N c . The mirror diagrams where the gluon attaches to the right of the cut 
line are not shown, but also contribute. 



consequence of the final state interaction effects in these SSAs. 

The rest of the paper is organized as follows. In Sec. 2, we will set up our framework, and 
calculate the SSAs for heavy quark and antiquark in the twist-three quark-gluon correlation 
approach and the initial/final state interactions effects are properly taken into account. 

We further present some numerical results for experiments at low energy pp scattering where 
the quark- antiquark channel dominates, and we will observe a factor of 3 difference between 
heavy quark and antiquark SSAs. We conclude our paper in Sec. 3. 

2. SSAs in Heavy Quark and Antiquark Productions. In order to obtain a non- 
zero SSA in high energy process, we must generate a phase from parton re-scattering in 
hard process [2, y|. Similar to other hadronic processes, in the heavy quark production in 
pp collisions 2|, [14( , the partons involved in the re-scattering can come from the initial state 
or final state, and they are referred as the initial and final state interactions, respectively. 
Because of the different color charges they carried in the scattering process, the heavy quark 
and antiquark will have different effects from the final state interactions. The initial state 
interactions, on the other hand, are the same for both particles. In this paper, we will take 
the quark- antiquark annihilation channel as an example to demonstrate the unique conse- 
quence from the different final state interaction effects on the heavy quark and antiquark 
productions. The light quark production in this channel has been studied in . We will 

extend to the heavy quark production and study the relevant phenomenological consequence. 

In the quark-antiquark annihilation channel, the quark (or antiquark) from the polarized 
nucleon annihilates the antiquark (or quark) from the unpolarized nucleon, and produces a 
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heavy quark and antiquark pair. For the spin-average cross section, this forms the s— channel 
partonic diagram. As we mentioned above, the single transverse spin dependent cross sec- 
tion comes from both initial and final state interaction contributions. In Fig. 1 we show the 
relevant initial (a) and final state interactions (b,c) diagrams. Different from those in the 
SIDIS and Drell-Yan lepton pair production, these initial and final state interactions have 
different strength because they have different color charges. This will result into different 
color-factors for them, which can be carried out in a model-independent way. For example, 
for the initial state interaction of Fig. 1(a), we can label the incoming quark legs from the 
polarized nucleon with color-index % and j ((i, j) = 1,3 in the elementary representation in 
the SU(3) gauge group), the additional gluon with index a (a = 1,8 for the adjoint repre- 
sentation) connecting the nucleon remanet with the partons in the hard partonic processes. 
The color-matrix representing the coupling of the nucleon with the partonic part will be 
uniquely the SU(3) color matriices T?j, because of the gauge invariance. Therefore, the 
color-factor associated with the initial state interaction can be formulated as 

1 M 2 — 1 

Fig.l(a) oc (-)tr [ T a T b T a T c ] x tr [T b T c ] = — c - , (1) 

where the additional minus sign is due to the initial state interaction effect. Similarly, the 
final state interaction diagram of Fig. 1(b) will be 

AT 2 — 9 /V 2 — 1 

Fig. 1(b) oc tr [T a T b T c ] x tr [T b T a T c ] = — . ( 2 ) 

Because this final state interaction is on the quark line, there is no additional minus sign. 
However, the final state interaction on the antiquark line will have a minus sign, 

2 /V 2 — 1 

Fig.l(c) oc (-)tr [T a T b T c ] x tr [ T a T b T c ] = c - . (3) 

From the above results, we find that all the initial and final state interactions contribute the 
same sign. However, the strengths of these interactions are different. Especially, the two 
final state interactions differ by a factor of 3. 

Furthermore, the initial state interaction diagram Fig. 1(a) contributes to the SSAs in 
both heavy quark and antiquark productions. However, the final state interaction on the 
quark line Fig. 1(b) only contributes to the heavy quark SSA, whereas that on the antiquark 
line Fig. 1(c) only contributes to that for the antiquark. For example, when we study the 
SSA for heavy quark production, the kinematics of heavy antiquark is integrated out. The 
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phase contribution from the final state interaction on the antiquark line (Fig. 1(c)) will cancel 
out that from the same diagram but with the gluon attachment to the right of the cut line 
(the mirror diagram of Fig. 1(c)) ^j. Therefore, the heavy quark SSA has contribution from 
the initial state interaction of Fig. 1(a) and the final state interaction on the quark line 
Fig. 1(b), whereas the heavy antiquark SSA from the initial state interaction of Fig. 1(a) 
and the final state interaction on the antiquark line Fig. 1(c). If the final state interactions 
dominate their SSAs, we will expect a factor of 3 difference between the SSAs in heavy 
quark and antiquark productions, because of the different strengths they have as we have 
shown in the above analysis. This is a model independent observation. Its test shall provide 
a strong support for the underlying physics for the SSA phenomena. 

To calculate the SSA contributions from these diagrams, it is appropriate to adopt the 
twist-three quark-gluon-antiquark correlation approach developed in [sj, because the heavy 
quark mass provides a natural hard scale in the computation of the perturbative diagrams, 
and the collinear factorization approach is applicable. An important feature from this ap- 
proach is that the single spin asymmetry is naturally suppressed in the limit of P± <C Mq, 
where P± and Mq are transverse momentum and mass of the heavy quark. 

nn 

The calculations are similar to those in [2|, H4|], where the inclusive it production in p'p — > 
irX was formulated. The only difference is that we have to restore the mass dependence for 
the final state quark and antiquark. The spin-average cross section can be summarized as, 

E ^ Q = t I "M f -ir *s + F-K) HS ^» l ' > l ' fi) • (4) 

where Pq is the heavy quark momentum, S the hadronic center of mass energy, S = (P4 + 
Pb) 2 with Pa and Pb the momenta for the incident polarized and unpolarized nucleons, x 
and x' are the momentum fractions carried by the incident partons, and f a and /& are the 
associated parton distributions. The kinematic variables are defined as: En = cosh y 



where Mt = y Mq + ^jl anc ^ V ls ^ ne ra pidity for the heavy quark in the center of mass 
frame; T = M 2 Q - M T y/Se~ y and U = M J - M T \fSe v ; the partonic variable s = xx'S, t = 
-xM T VSe- y , and u = —x' MT\/~Se y . The leading contribution to heavy quark production 
comes from quark-antiquark and gluon-gluon channels, and their hard factor H u are defined 
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clS, 

v J_/J_ 2iV c 2 1 \ feg + n 2 + 4M 2 ) - 4g 2 M* 

^-QQ 4iV c ^~n TV 2 -Is 2 J tu ' [) 

In the twist-three framework Q], the single transverse spin dependent cross section depends 
on the twist-three quark-gluon correlation function, the so-called Qiu-Sterman matrix ele- 
ment. They are defined as 



d( dr]~ 

47T 

PS\%(o)£(ox-h + gF + a (nc(C-,v-)Mir)\PS) , 



n( Xl ,x 2 ) = / ^L e ^-^C-) e ^ S±p (7) 



where the sums over color and spin indices are implicit, e? a = e^ a P AfjL P Bu / P A ■ P B with 
e 0123 = 1, \PS) denotes the proton state, ip the quark field, F + a the gluon field tensor, and 
£ represents the gauge link to make the above definition gauge invariant. By summing the 
initial and final state interaction contributions, the single spin dependent differential cross 
section for heavy quark production can be written as 

d*Aa(S ± ) of /-dx' e^SlPj 1 f (T F (x) d 

d*P Q S J x' 5 x'S + (T - M^) u\\ x dx F[ ' 

T q (x) ~ 1 

X H qq ^Q + — H q q_Q > , (8) 



x 

where T F (x) = T F (x, x). A similar expression can be obtained for the antiquark by replacing 
Q — > Q in the above. The above hard factors include both initial and final state interaction 
contributions. The first term H qq -_>Q follows a compact formula containing the derivative 



and non-derivative (of the twist-three function) contributions [14|, whereas the second term 
H q q_,Q only depends on the non-derivative of the twist-three correlation function and van- 
ishes in the massless limit Mn — * 0. Under this limit, our results reproduce the relevant 

n 

formula in |14J]. The hard factors are found, 

H qq ^Q = Hq q ^ Q + Hqq^ Q (l + , (9) 

u , 



Hqq^Q = Hqq^Q + H^ Q (1 + j) , (10) 
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where 



1 2(t 2 + u 2 + 2Mq§) ~ r ~ r I Ml 

ttI _ ttI _ _ v \ Q 1 ttI _ ttI _ _ x Q_ 

QQ—+Q ~ n aa^O~ A at-?. ~9 ' n <W-><2 ~~ n aa-^O ~ at-?. ~ ' V 11 / 



qq^Q qq-*Q 4^2 £2 ' qq^Q qq-^Q jy2 g 

2 2 (t 2 + « 2 + 2M 2 s) 2 M 2 



iV 2 -2 2(P + S 2 + 2M 2 S) - F _iV 2 -2M 2 

n QQ^O. ~ A AT9. ~9 ' n qq^Q. ~ AT9. ~ ' ^ 1D J 



-qq-^Q AN 2 ~ s 2 ' qq-MJ ft* g 

where the color factors have been discussed in the above. The contributions from the 
antiquark-gluon correlation functions can be obtained accordingly, by using the charge con- 
jugation transformation. 

It is important to note that the spin-average and spin-dependent cross section contri- 
butions are both well defined in the limit of P± — > 0, because of the heavy quark mass. 
Therefore, the spin asymmetry defined as the ratio of these two cross section terms vanishes 
when P± = 0. This is very different from the massless particle production where the single 
spin asymmetry divergent as 1/P± as P± — > . Thus we can integrate out the cross 

sections to small transverse momentum region to obtain the spin asymmetry. In the follow- 
ing, we will present the numerical results for the conventional left-right asymmetry for the 
heavy quark production at low energy pp scattering. This asymmetry is defined as 

(14) 

which counts the number difference between the left hand side and right hand side for heavy 
quark production from the moving direction of the polarized (upward) nucleon. It can be 
calculated from the cross section terms from above, 

An ~ £d<pda ' (15) 

where <ft is the azimuthal angle of heavy quark transverse momentum P± relative to the 
incoming nucleon transverse polarization vector S±. In the following, we will show the 
above asymmetry for heavy quark and antiquark at low energy pp scattering where the 
quark- ant iquark channel dominates. For the purpose of the demonstration, we adopt the 
paramterizations for the twist-three Qiu-Sterman matrix elements in the valence region 
as Q, 

T a F (x) = N a x a «(l-xf«q a (x) , (16) 
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FIG. 2: Single spin asymmetries for charm quark and anticharm quark production inp^p scattering 
at JPARC energy: as functions of rapidities with P± integrated out (left panel); and as functions 
of transverse momentum P± with y c .m. > (right panel). 

where q a (x) is the unpolarized quark distribution for flavor a. The above parameters for the 
valence u and d quarks were determined from a comparison to the single spin asymmetry data 



in p^p 
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141 ]. and they are also compatible with the SIDIS data from HERMES 15] 



For the unpolarized parton distributions, we use the leading order (LO) CTEQ5L parton 
distribution functions [l^ ]. 

In Fig. 2, we show our predictions for the heavy quark and antiquark single spin asym- 
metries in p^p scattering at JPARC energy region (y/s ~ 7.8GeV), where we expect the 
quark- antiquark annihilation channel dominates the cross section. The left panel shows the 
asymmetries as functions of the center of mass rapidities of the charm quark and anti- 
charm quark with the transverse momentum P± integrated out; the right panel shows the 
asymmetries as functions of the transverse momentum P± in the forward rapidity region 
(Ucm. > 0). From these plots, we find that the asy mmetries increase with rapidity, similar to 
the general trend in other SSA observations 14j. More importantly, the single spin asym- 
metry for the charm quark is much larger than that for anticharm quark in the full rapidity 
range. Similar observation was found for the asymmetries as functions of the transverse 
momentum P±. Both asymmetries vanish as Pj_ — ► as we expected. 

In Fig. 3, we show the predictions at \/s = 14GeV p^p scattering at the PAX experiment 
at GSI PAIR. Because the quark-antiquark channel always dominates the cross section, we 
find a sizable SSA for charm quark even at cental rapidity region. Again, we observe that 
the SSA for charm quark is about a factor of 3 larger than that for the anticharm quark. In 
this experiment, we can further study the spin asymmetry with the antiproton transversely 
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FIG. 3: Single spin asymmetries for charm and anticharm quark production in p^p scattering at 
y/s = \AGeV at PAX energy region: as functions of raptidities with P± integrated out (left panel); 
as functions of charm quark transverse momentum P± with y c . m . > (right panel). 





FIG. 4: Ratio of the SSA for the charm quark over that for the anticharm quark as function of charm 
quark rapidity inp^p scattering at JPARC (left) andp^p scattering at PAX (right) experiments. 

polarized p^p — > Q(Q)X, where we expect the SSA for anticharm quark is larger than that 
for the charm quark in the forward direction of the polarized antiproton. By using the 
charge conjugation transformation invariance, we will have the following relation between 
the SSAs in these two processes, 

A N tfp - C) = A N (plp - C) , A N tfp - C) = A N (plp - C) , (17) 

for the same kinematics. From these relations, we can obtain the corresponding results for 
p^p experiments from the plots in Fig. 3. 

The sizes of the above single spin asymmetries also depend on the dominance of the 
quark- ant iquark annihilation channel at both energies. In our leading order simulation with 
the CTEQ5L parameterization for the parton distributions, this channel indeed dominates 
the cross sections. These asymmetries may scale down if the gluon-gluon channel turns 
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dominant 



17| . However, we shall still observe the big difference between the SSAs for 



charm and anticharm quarks when the twist-three gluon-gluon correlation contribution to 
the SSAs 3, is small in the valence region. In order to demonstrate this difference, we 
plot in Fig. 4 the ratios of the asymmetries as functions of the rapidities at both experiments. 
From these plots, we find that the ratio at both experiments is about a factor of 3 and has 
little dependence on the rapidity. These ratios are model-independent predictions, especially 
at forward region of the polarized nucleon where the quark-antiquark channel dominates. 
We hope the experiments at both facilities can carry out these measurements, and test the 
predictions. We further notice that the sea quark contribution in the quark-gluon correlation 
functions from the polarized nucleon will be important at the central rapidities [14J], where 
the above ratios may vary if we take into account those contributions. 

4. Summary. In this paper, we studied the single spin asymmetries in heavy quark 
and antiquark production in hadronic processes. Because of the different color charges, the 
final state interactions on the heavy quark and antiquark are different, and therefore the 
associated single spin asymmetries are different too. Our numerical simulations showed that 
the single spin asymmetry for heavy quark is about a factor of 3 larger than that for the 
antiquark in the forward p^p amd p^p scattering at low energy sacttering where the quark- 
antiquark channel dominates. The experimental observation of this unique signature shall 
provide a crucial test for the underlying physics for the single spin asymmetry phenomena. 

In the above analysis, we only considered the quark-antiquark channel contributions. An 
extension to the gluon-channel contribution can follow accordingly, which will be relevant for 
the collider experiment at RHIC 18, Q, [2^]. We reserve the study of this contribution as well 



as other contribution, for example, from the chiral-odd quark-gluon correlation functions 
in a future publication. We further notice that the latter contribution will have the similar 
pattern for the SSAs in heavy quark and antiquark production as the chiral-even one as 
we discussed in this paper, because of the final state interaction effects dominate both 
contributions. 
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